The scale of population structure in Arabidopsis thaliana.

The population structure of an organism reflects its evolutionary history and influences its evolutionary trajectory. It constrains the combination of genetic diversity and reveals patterns of past gene flow. Understanding it is a prerequisite for detecting genomic regions under selection, predicting the effect of population disturbances, or modeling gene flow. This paper examines the detailed global population structure of Arabidopsis thaliana. Using a set of 5,707 plants collected from around the globe and genotyped at 149 SNPs, we show that while A. thaliana as a species self-fertilizes 97% of the time, there is considerable variation among local groups. This level of outcrossing greatly limits observed heterozygosity but is sufficient to generate considerable local haplotypic diversity. We also find that in its native Eurasian range A. thaliana exhibits continuous isolation by distance at every geographic scale without natural breaks corresponding to classical notions of populations. By contrast, in North America, where it exists as an exotic species, A. thaliana exhibits little or no population structure at a continental scale but local isolation by distance that extends hundreds of km. This suggests a pattern for the development of isolation by distance that can establish itself shortly after an organism fills a new habitat range. It also raises questions about the general applicability of many standard population genetics models. Any model based on discrete clusters of interchangeable individuals will be an uneasy fit to organisms like A. thaliana which exhibit continuous isolation by distance on many scales.

Salicylic acid biosynthesis is enhanced and contributes to increased biotrophic pathogen resistance in Arabidopsis hybrids.

Heterosis, the phenotypic superiority of a hybrid over its parents, has been demonstrated for many traits in Arabidopsis thaliana, but its effect on defence remains largely unexplored. Here, we show that hybrids between some A. thaliana accessions show increased resistance to the biotrophic bacterial pathogen Pseudomonas syringae pv. tomato (Pst) DC3000. Comparisons of transcriptomes between these hybrids and their parents after inoculation reveal that several key salicylic acid (SA) biosynthesis genes are significantly upregulated in hybrids. Moreover, SA levels are higher in hybrids than in either parent. Increased resistance to Pst DC3000 is significantly compromised in hybrids of pad4 mutants in which the SA biosynthesis pathway is blocked. Finally, increased histone H3 acetylation of key SA biosynthesis genes correlates with their upregulation in infected hybrids. Our data demonstrate that enhanced activation of SA biosynthesis in A. thaliana hybrids may contribute to their increased resistance to a biotrophic bacterial pathogen.

Cell Fate Specification and the Regulation of RNA-dependent DNA Methylation in the Arabidopsis Root Meristem

The Arabidopsis root apical meristem (RAM) is a complex tissue capable of generating all the cell types that ultimately make up the root. The work presented in this thesis takes advantage of the versatility of high-throughput sequencing to address two independent questions about the root meristem. Although a lot of information is known regarding the cell fate decisions that occur at the RAM, cortex specification and differentiation remain poorly understood. In the first part of this thesis, I used an ethylmethanesulfonate (EMS) mutagenized marker line to perform a forward genetics screen. The goal of this screen was to identify novel genes involved in the specification and differentiation of the cortex tissue. Mapping analysis from the results obtained in this screen revealed a new allele of BRASSINOSTEROID4 with abnormal marker expression in the cortex tissue. Although this allele proved to be non-cortex specific, this project highlights new technology that allows mapping of EMS-generated mutations without the need to map-cross or back-cross. In the second part of this thesis, using fluorescence activated cell sorting (FACS) coupled with high throughput sequencing, my collaborators and I generated single-base resolution whole genome DNA methylomes, mRNA transcriptomes, and smallRNA transcriptomes for six different populations of cell types in the Arabidopsis root meristem. We were able to discover that the columella is hypermethylated in the CHH context within transposable elements. This hypermethylation is accompanied by upregulation of the RNA-dependent DNA methylation pathway (RdDM), including higher levels of 24-nt silencing RNAs (siRNAs). In summary, our studies demonstrate the versatility of high-throughput sequencing as a method for identifying single mutations or to perform complex comparative genomic analyses.

Information processing without brains--the power of intercellular regulators in plants.

Plants exhibit different developmental strategies than animals; these are characterized by a tight linkage between environmental conditions and development. As plants have neither specialized sensory organs nor a nervous system, intercellular regulators are essential for their development. Recently, major advances have been made in understanding how intercellular regulation is achieved in plants on a molecular level. Plants use a variety of molecules for intercellular regulation: hormones are used as systemic signals that are interpreted at the individual-cell level; receptor peptide-ligand systems regulate local homeostasis; moving transcriptional regulators act in a switch-like manner over small and large distances. Together, these mechanisms coherently coordinate developmental decisions with resource allocation and growth.

Detecting separate time scales in genetic expression data.

BACKGROUND: Biological processes occur on a vast range of time scales, and many of them occur concurrently. As a result, system-wide measurements of gene expression have the potential to capture many of these processes simultaneously. The challenge however, is to separate these processes and time scales in the data. In many cases the number of processes and their time scales is unknown. This issue is particularly relevant to developmental biologists, who are interested in processes such as growth, segmentation and differentiation, which can all take place simultaneously, but on different time scales. RESULTS: We introduce a flexible and statistically rigorous method for detecting different time scales in time-series gene expression data, by identifying expression patterns that are temporally shifted between replicate datasets. We apply our approach to a Saccharomyces cerevisiae cell-cycle dataset and an Arabidopsis thaliana root developmental dataset. In both datasets our method successfully detects processes operating on several different time scales. Furthermore we show that many of these time scales can be associated with particular biological functions. CONCLUSIONS: The spatiotemporal modules identified by our method suggest the presence of multiple biological processes, acting at distinct time scales in both the Arabidopsis root and yeast. Using similar large-scale expression datasets, the identification of biological processes acting at multiple time scales in many organisms is now possible.

Intergenic and genic sequence lengths have opposite relationships with respect to gene expression.

Eukaryotic genomes are mostly composed of noncoding DNA whose role is still poorly understood. Studies in several organisms have shown correlations between the length of the intergenic and genic sequences of a gene and the expression of its corresponding mRNA transcript. Some studies have found a positive relationship between intergenic sequence length and expression diversity between tissues, and concluded that genes under greater regulatory control require more regulatory information in their intergenic sequences. Other reports found a negative relationship between expression level and gene length and the interpretation was that there is selection pressure for highly expressed genes to remain small. However, a correlation between gene sequence length and expression diversity, opposite to that observed for intergenic sequences, has also been reported, and to date there is no testable explanation for this observation. To shed light on these varied and sometimes conflicting results, we performed a thorough study of the relationships between sequence length and gene expression using cell-type (tissue) specific microarray data in Arabidopsis thaliana. We measured median gene expression across tissues (expression level), expression variability between tissues (expression pattern uniformity), and expression variability between replicates (expression noise). We found that intergenic (upstream and downstream) and genic (coding and noncoding) sequences have generally opposite relationships with respect to expression, whether it is tissue variability, median, or expression noise. To explain these results we propose a model, in which the lengths of the intergenic and genic sequences have opposite effects on the ability of the transcribed region of the gene to be epigenetically regulated for differential expression. These findings could shed light on the role and influence of noncoding sequences on gene expression.

The Function and Regulation of Photobodies in Phytochrome Signaling

Light is a critical environmental signal that regulates every phase of the plant life cycle, from germination to floral initiation. Of the many light receptors in the model plant Arabidopsis thaliana, the red- and far-red light-sensing phytochromes (phys) are arguably the best studied, but the earliest events in the phy signaling pathway remain poorly understood. One of the earliest phy signaling events is the translocation of photoactivated phys from the cytoplasm to the nucleus, where they localize to subnuclear foci termed photobodies; in continuous light, photobody localization correlates closely with the light-dependent inhibition of embryonic stem growth. Despite a growing body of evidence supporting the biological significance of photobodies in light signaling, photobodies have also been shown to be dispensable for seedling growth inhibition in continuous light, so their physiological importance remains controversial; additionally, the molecular components that are required for phy localization to photobodies are largely unknown. The overall goal of my dissertation research was to gain insight into the early steps of phy signaling by further defining the role of photobodies in this process and identifying additional intragenic and extragenic requirements for phy localization to photobodies.

Even though the domain structure of phys has been extensively studied, not all of the intramolecular requirements for phy localization to photobodies are known. Previous studies have shown that the entire C-terminus of phys is both necessary and sufficient for their localization to photobodies. However, the importance of the individual subdomains of the C-terminus is still unclear. For example a truncation lacking part of the most C-terminal domain, the histidine kinase-related domain (HKRD), can still localize to small photobodies in the light and behaves like a weak allele. However, a point mutation within the HKRD renders the entire molecule completely inactive. To resolve this discrepancy, I explored the hypothesis that this point mutation might impair the dimerization of the HKRD; dimerization has been shown to occur via the C-terminus of phy and is required for more efficient signaling. I show that this point mutation impairs nuclear localization of phy as well as its subnuclear localization to photobodies. Additionally, yeast-two-hybrid analysis shows that the wild-type HKRD can homodimerize but that the HKRD containing the point mutation fails to dimerize with both itself and with wild-type HKRD. These results demonstrate that dimerization of the HKRD is required for both nuclear and photobody localization of phy.

Studies of seedlings grown in diurnal conditions show that photoactivated phy can persist into darkness to repress seedling growth; a seedling's growth rate is therefore fastest at the end of the night. To test the idea that photobodies could be involved in regulating seedling growth in the dark, I compared the growth of two transgenic Arabidopsis lines, one in which phy can localize to photobodies (PBG), and one in which it cannot (NGB). Despite these differences in photobody morphology, both lines are capable of transducing light signals and inhibiting seedling growth in continuous light. After the transition from red light to darkness, the PBG line was able to repress seedling growth, as well as the accumulation of the growth-promoting, light-labile transcription factor PHYTOCHROME INTERACTING FACTOR 3 (PIF3), for eighteen hours, and this correlated perfectly with the presence of photobodies. Reducing the amount of active phy by either reducing the light intensity or adding a phy-inactivating far-red pulse prior to darkness led to faster accumulation of PIF3 and earlier seedling growth. In contrast, the NGB line accumulated PIF3 even in the light, and seedling growth was only repressed for six hours; this behavior was similar in NGB regardless of the light treatment. These results suggest that photobodies are required for the degradation of PIF3 and for the prolonged stabilization of active phy in darkness. They also support the hypothesis that photobody localization of phys could serve as an instructive cue during the light-to-dark transition, thereby fine-tuning light-dependent responses in darkness.

In addition to determining an intragenic requirement for photobody localization and further exploring the significance of photobodies in phy signaling, I wanted to identify extragenic regulators of photobody localization. A recent study identified one such factor, HEMERA (HMR); hmr mutants do not form large photobodies, and they are tall and albino in the light. To identify other components in the HMR-mediated branch of the phy signaling pathway, I performed a forward genetic screen for suppressors of a weak hmr allele. Surprisingly, the first three mutants isolated from the screen were alleles of the same novel gene, SON OF HEMERA (SOH). The soh mutations rescue all of the phenotypes associated with the weak hmr allele, and they do so in an allele-specific manner, suggesting a direct interaction between SOH and HMR. Null soh alleles, which were isolated in an independent, tall, albino screen, are defective in photobody localization, demonstrating that SOH is an extragenic regulator of phy localization to photobodies that works in the same genetic pathway as HMR.

In this work, I show that dimerization of the HKRD is required for both the nuclear and photobody localization of phy. I also demonstrate a tight correlation between photobody localization and PIF3 degradation, further establishing the significance of photobodies in phy signaling. Finally, I identify a novel gene, SON OF HEMERA, whose product is necessary for phy localization to photobodies in the light, thereby isolating a new extragenic determinant of photobody localization. These results are among the first to focus exclusively on one of the earliest cellular responses to light - photobody localization of phys - and they promise to open up new avenues into the study of a poorly understood facet of the phy signaling pathway.

Redox rhythm reinforces the circadian clock to gate immune response.

Recent studies have shown that in addition to the transcriptional circadian clock, many organisms, including Arabidopsis, have a circadian redox rhythm driven by the organism's metabolic activities. It has been hypothesized that the redox rhythm is linked to the circadian clock, but the mechanism and the biological significance of this link have only begun to be investigated. Here we report that the master immune regulator NPR1 (non-expressor of pathogenesis-related gene 1) of Arabidopsis is a sensor of the plant's redox state and regulates transcription of core circadian clock genes even in the absence of pathogen challenge. Surprisingly, acute perturbation in the redox status triggered by the immune signal salicylic acid does not compromise the circadian clock but rather leads to its reinforcement. Mathematical modelling and subsequent experiments show that NPR1 reinforces the circadian clock without changing the period by regulating both the morning and the evening clock genes. This balanced network architecture helps plants gate their immune responses towards the morning and minimize costs on growth at night. Our study demonstrates how a sensitive redox rhythm interacts with a robust circadian clock to ensure proper responsiveness to environmental stimuli without compromising fitness of the organism.

The origin and evolution of phototropins.

Plant phototropism, the ability to bend toward or away from light, is predominantly controlled by blue-light photoreceptors, the phototropins. Although phototropins have been well-characterized in Arabidopsis thaliana, their evolutionary history is largely unknown. In this study, we complete an in-depth survey of phototropin homologs across land plants and algae using newly available transcriptomic and genomic data. We show that phototropins originated in an ancestor of Viridiplantae (land plants + green algae). Phototropins repeatedly underwent independent duplications in most major land-plant lineages (mosses, lycophytes, ferns, and seed plants), but remained single-copy genes in liverworts and hornworts-an evolutionary pattern shared with another family of photoreceptors, the phytochromes. Following each major duplication event, the phototropins differentiated in parallel, resulting in two specialized, yet partially overlapping, functional forms that primarily mediate either low- or high-light responses. Our detailed phylogeny enables us to not only uncover new phototropin lineages, but also link our understanding of phototropin function in Arabidopsis with what is known in Adiantum and Physcomitrella (the major model organisms outside of flowering plants). We propose that the convergent functional divergences of phototropin paralogs likely contributed to the success of plants through time in adapting to habitats with diverse and heterogeneous light conditions.